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Moduli evolution in the presence of flux compactifications 
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We study the cosmological evolution of the volume moduli in a class of recently proposed Infla- 
tionary Universe models arising out of Type IIB string theory where a number of the moduli 
fields have been stabilised through flux compactifications 0. Developing an approach initially in- 
troduced in Q we show, in agreement with Q, how the presence of extra sources of matter act so 
as to provide additional friction, slowing the modulus field as it evolves down its potential, thereby 
vastly increasing the region of parameter space which leads to the eventual stabilisation of these 
fields. Extending the case to include both the real and imaginary parts of the volume modulus, we 
show how the parameter space of inital conditions is modified and comment on the impact for these 
inflationary models arising out of flux type compactifications. 



I. INTRODUCTION 

String theory has associated with it many moduli 
fields, scalar fields whose presence are important in deter- 
mining quantities such as the size of the internal dimen- 
sions, the gauge coupling constants, even the strength 
of the gravitational interactions. One of the problems 
that has faced physicists attempting to develop cosmol- 
ogy within string theory has been how to deal with these 
fields - basically they never want to stop evolving. In 
most models that have been developed to date, the mod- 
uli initially have no associated potential with which to 
trap them, the flat directions mean that there is noth- 
ing to prevent them from evolving for ever, leading to 
time varying coupling constants, decompactification of 
the internal dimensions and other nightmares. Non triv- 
ial potentials have been shown to emerge when non- 
perturbative features are included, such as gaugino con- 
densates but even then the local minima associated 
with these models are generally unstable in that the bar- 
rier height protecting it from the runaway solution is 
small compared to any other natural scale in the prob- 
lem. Moreover, the value of the potential at its minima 
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is negative, implying an Anti de Sitter Space solution |g] . 
This issue has plagued string models for many years. Ini- 
tially, in the context of weakly coupled heterotic string 
theory, this was pointed out by Dine and Seiberg 0. In 
the context of cosmology Brustein and Steinhardt read- 
dressed the issue, pointing out how difficult it would be 
to stabilise typical moduli fields emerging in these models 
because of the combination of the steepness of the poten- 
tials associated with the moduli, and the small size of the 
potential barrier separating the fields from the runaway 
regime §|. Earlier, in ref. 0, Kaloper and Olive had 
studied the effect including a background of radiation in 
the equations that govern the dilaton evolution. In par- 
ticular, they found that the dilaton can stop rolling if 
cither its potential has a SUSY breaking minimum or ra- 
diation is dominating the energy density of the Universe. 
In Q this idea was further developed, in particular it 
was pointed out that, in the presence of an extra fluid 
component dominating the energy density, for example 
radiation, the effect of the fluid was to increase the effec- 
tive friction experienced by the moduli field as it evolved 
down its steep potential. For a wide range of initial con- 
ditions of the field it was demonstrated that it would be 
slowed down, and its energy density would tend towards 
a scaling solution where it followed the evolution of the 
background radiation energy density. Moreover, it would 
then simply fall into the minimum of the potential and 
be stabilised there, even though the height of the poten- 
tial barrier was so low. This analysis, which was initially 
used to demonstrate how the dilaton could be stabilised 
in weakly coupled heterotic string theory, was later ex- 
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tended to the case of heterotic M-theory, where including 
the evolution of more than one modulus field, it became 
evident how important the coupling between the various 
moduli fields could be in determining the basin of attrac- 
tion into the true minimum of the potential [10( . Other 
attempts at stabilising the moduli fields in these models 
included adding specific temperature corrections to the 
moduli fields, in order to extend the range of initial val- 
ues of the fields which would lead to stabilisation [Tlj . 
although this met with limited success [ljj . 

Over the past few years there has been a renewed in- 
terest in the issue of moduli evolution. Giddings et al 
demonstrated, in the context of Type IIB string theory 
compactified on a Calabi-Yau manifold, how it is possible 
to use fluxes to stabilise all but the Kahler moduli, in par- 
ticular the overall volume modulus [l^. This result was 
used to develop inflation models in string theory where 
brane inflation was embedded into one of these stable 
compactified models 0, Q ■ However in pj it was realised 
that there was a price to pay, and a degree of fine tuning 
was required in these models in order to obtain sufficient 
inflation. The particular problem they encountered was 
that, as the compactification volume modulus was being 
stabilised, the effect was to modify the inflaton potential, 
rendering it too steep for inflation to last for a sufficient 
length of time. The issue of stabilisation had raised its 
ugly head again. Recently Brustcin ct al Q returned to 
the idea of using sources as a way of increasing the fric- 
tion and thereby avoiding the runaway problem alluded 
to earlier. They did it in the context of a toy model with 
a single scalar field (the volume modulus) and scalar po- 
tential of the form suggested in [|. They argued that 
the key reason they could stabilise the modulus was be- 
cause of the way this field quickly entered a regime where 
its energy density was dominated by its kinetic energy 
It therefore lost its energy far quicker than any other 
source, allowing the other sources (such as radiation) to 
eventually catch up and take over, slowing the field down 
even more and allowing it to come to rest at or near the 
local minimum of its potential. In this paper we extend 
the investigation into these models by allowing for the 
fact that the volume modulus is a complex field with 
a real and imaginary component. We therefore look at 
the system where more than one moduli field has to be 
dynamically stabilised in the presence of extra sources. 
We reproduce the results of 3 in the appropriate limits, 
but show that there are wider ranges of initial conditions 
which can lead to stability if we include the possibility 
of the moduli fields entering a period of scaling, where 
their energy density mimics that of the background en- 
ergy density, more in the spirit of @ E3- A related 
approach to the stabilisation issue has recently been in- 
vestigated in refs. 0,0,0,0, where the authors have 
considered the presence of a gas of wrapped branes on the 
dynamics of some of the moduli fields. 

The layout of this paper is as follows. In section ITT1 we 
introduce the system of equations we will be investigat- 
ing in terms of a general moduli potential, and briefly 



describe the nature of the attractor solutions when ex- 
ponential terms are present in the scalar potential. The 
particular example due to Kachru et al p| is analysed 



in section UTll and a related but different example due to 
Kallosh and Linde 0] is examined in section lTvl Finally 
we conclude in section IV! 



II. BACKGROUND EQUATIONS OF MOTION 

The class of models that we will be investigating can 
in general be described by an N=l, d=4 effective Su- 
pergravity (SUGRA) theory, where d is the number of 
non-compact spacetime dimensions. In this case the four- 
dimensional N=l SUGRA action is of the form 
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. is the Kahler metric: $ 4 are complex 
chiral supcrficlds; V(<f>) is the scalar potential and Kp is 
the 4-dimensional Newton constant which, in terms of 
higher dimensional quantities, can be expressed as 



2irpv 



N 



(2) 



where n is the 11-dimcnsional Newton constant, v is the 
volume of the six dimensional compact manifold, and p 
is the radius of the eleventh spatial dimension. 

The impact of the non-perturbative and flux effects 
considered in Q, is to induce an effective scalar poten- 
tial for the volume moduli whose most general form for 
4-dimensional N = 1 SUGRA is 



V = e K (K i: >D l WD j W - 3WW) 



(3) 



where K % $ is the inverse Kahler metric and DiW = 
diW + §§tW is the Kahler covariant derivative acting 
on the superpotential. 

The equations of motion follow from the variation of 
the action (JJ. For simplicity we consider the case of 
homogeneous, timc-dependent, fields in a spatially flat 
Fricdmann-Robertson- Walker space time background. 
Given this ansatz, we find the following equations of mo- 
tion for the complex superfields 



(4) 



$' + 3H& + r ! ; ,,$ 3 $ fc + K %3 d-Y = , 

where are the relevant moduli, $ 4 = d^/dt, djV = 

dV/d§ 3 , and the connection on the Kahler manifold has 
the form 



pn jy-nt J 1 



(5) 



In addition, we obtain the Friedman equation for the 
Hubble factor H = a/a, where a(t) is the scale factor of 
the Universe, 

3H 2 = 4(p 4 + Pb ) = K %(K i3 && + V + Pb ) , (6) 
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with p$ and pb are the energy density of the evolving 
moduli fields and background fluid respectively. The dy- 
namics of the latter is given in terms of the scale factor 
and its background equation of state, 7 — 1 = Pb/ Pb, 
where pb is the pressure of the fluid, 



Pb = Pbo/a 



3 7 



(7) 



In what follows we set n 2 p 



1. 



It is worth splitting the equations of motion for the 
complex chiral superfields into those for their real and 
imaginary parts 



( p R + 3Hf R +T) k (^ R -^ k I ) + ^d jR V = , (8) 
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4>\ + 3H4>} + r jk io-o'j, + 4P R 4>1) + ^ d n v = o , (9) 



where now 



refers to the real (imaginary) part 



of the scalar fields and dj R {dj, ) are used to denote the 
derivative of the potential with respect to the real (imag- 
inary) parts of the fields respectively 



The system of differential equations that one has to 
solve comes from eqs. <|8I9[) being applied to the present 
model, along with the evolution equation for the back- 
ground fluid. Altogether, we have 
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cfi + 3HJi d r di + —^d^V = , (13) 

oy 3 

p b + iH-fpb = , 
subject to the Friedman constraint, see eq. 10, 
3 



'iH 2 = —{d 2 +d 2 ) + V + Pb 
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The scalar potential acquires a relatively simple form 
when written in terms of both real and imaginary parts 
of a, 
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III. KKLT MODEL 

The possibility of finding dc Sitter vacua in string the- 
ory with a stabilized volume modulus, er, was put forward 
in rcf. 0, and has been widely adopted in subsequent 
work. The key ingredient was to consider the combina- 
tion of non perturbative effects and an additional flux 
term in the superpotential 

W = W + A&- aa , (10) 

which, together with the usual Kahlcr potential 



K = -31n((7 + a) , 



(11) 



defines the F-part of the SUGRA potential, see eq. (J3J). It 
has been known for many years now that, in this context, 
it is possible to stabilize a, although giving rise to an Anti 
de Sitter (AdS) vacuum. As pointed out in rcf. ,2], if wc 
include contributions from cither anti-D3 or D7 brancs, 
an additional D-type term is generated, of the form 



V D = ^ 



(12) 



where we write a = a r + ia^. By suitably tuning the 
value of C one can move to a dS -or even Minkowski- 
vacuum. 

In this section we are interested in studying the cos- 
mological evolution of the field a as it rolls towards its 
minimum. Previous results addressing the same issue 
were published in ref . , where only the evolution of the 
real part of a, o>, was considered. Here we would like to 
extend this to study the behaviour of both a r and Ui. 



Given the above expression it is easy to see that the po- 
tential has an extremum in <Ji for aoi = tit:, with n an 
integer. Depending on the sign of Wocos(a<ji) this can 
be either a maximum or a minimum. For example, in 
figure n we show a contour plot of this scalar potential, 




a a./it 



FIG. 1: Contour plot of the scalar potential, given by eq. 1151 1. 
in the (<j>, a<Ji/ir) plane (solid lines). The values of the param- 
eters are taken from ref. Q, see text, and the dashed lines cor- 
respond to the supersymmetry-preserving conditions, F a = 0, 
where F a = K a W + W„. 

V, in the ((p,(Ti) plane, for the same values used in ref. 0], 
namely A = 1.0, a = 0.1, C = 3 x 10~ 26 , and W neg- 
ative (with cos(aCTi) = 1) and such that the minimum 
at (Ji — is supersymmetric (see below). Following that 
same reference we work with the canonically normalized 
field = y / 3/ 7 21n<7 r , instead of a r itself. 
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It is also worth mentioning the supersymmctric char- 
acter of the minima shown in this plot. We have included 
the two supersymmetry-preserving conditions, Re(F CT ) = 
Im(Fcr) = 0, and it can be clearly seen how both meet 
at the minima. Those correspond to 4> = 7.06 and even 
values of aai/n. 



A. One field evolution 



In figurc|31we present the initial values of the variables 
(</>,f2 r ) for which there is stabilization of the field at the 
minimum of its potential in the presence of radiation p_). 
i.e. 7 = 4/3. As expected, the smaller the initial back- 
ground fraction (parametrized by f2 r ) is, the smaller the 
region of allowed initial values of <fi which lead to a late 
time stabilisation of the field at the minimum. The solid 
line along Q r = 0.5 corresponds to the allowed region 
quoted in ref. 



In order to discuss the results of ref. |4| , we will first 
set the potential at a minimum in oi and solve only the 
first and third equations of the system l|13|) . This would 
correspond to taking the slice aoi = in figure 1, and 
considering only the evolution along 0. As for the initial 
conditions, we set oy = 0, and we choose values for cr-o 
and the fractional energy density in the background fluid, 
fl b = K 2 pPb0 /3Hl We can then calculate the value of pto, 
using the relation 



Pba = V(a r0 ) 



n,, 



(16) 



In this section, the results are shown in terms of the initial 
abundance flh which will allow us to compare them with 
those in ref. [J, where fl b = 0.5. 

In general one can identify up to five regions in the 
evolution of a scalar field with these types of potentials. 
In figurc|21we show a typical evolution going through the 
five regions although, of course, not all initial conditions 
will give rise to an evolution that will go through all of 
them. 

First, if the energy density of the background domi- 
nates, the field is effectively frozen in its evolution. This 
can be seen in the plot before region 1, with an evolution 
similar to region 3. When the energy density in the back- 
ground becomes comparable (region 1), the field starts 
rolling down its potential and eventually dominates the 
dynamics of the Universe. This happens when the po- 
tential is very shallow and is called the scalar field dom- 
inated solution. Eventually the potential becomes very 
steep, leading the evolution to a kinetically dominated 
solution (region 2), which ends when the frictional term 
in the equation of motion becomes dominant . The length 
of this period of "kination" is related to the value of the 
ratio Pb/P4> at the end of region 1. In region 3 the field 
is effectively frozen in the potential. The field restarts 
rolling down the potential once the background energy 
density has decayed to a value such that the frictional 
and potential terms in the equation of motion balance 
each other. At this stage the field evolves with a nearly 
constant ratio between kinetic and potential energy. This 
is called the scaling (or tracker) solution (region 4). Fi- 
nally, if this ratio is sufficiently small the field does not 
possess enough kinetic energy to roll over the potential 
barrier and gets traped at the minimum (region 5). De- 
tails of all these phases of the evolution are explained in 
Appendix A. 
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FIG. 2: Evolution of the energy densities of the field 4> (solid 
line) and the background (dashed line) with the logarithm of 
the scale factor. The various epochs of the evolution, num- 
bered from 1 to 5, are described in the text. We have taken 
0o = —5, 7 = 1 and Q m = 0.99. The parameters of the scalar 
potential are the same as used in ref. |1, namely A = 1.0, 
a = 0.1, C = 3 x 10" 26 . 
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FIG. 3: Fraction of the total initial energy needed in radia- 
tion, as a function of the initial value of <j>, for the field to end 
up in its minimum. The solid line along Q, r = 0.5 denotes the 
result obtained in ref. 3 and explained in the text. 
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B. Two field evolution 

Let us now introduce cr,, the imaginary part of the 
scalar field a in the evolution. The result of solving the 
full system given by eqs. (|13[) is presented in figure El 
where we show contour plots for the initial values of the 
pair (</>, a,), that lead to both fields ending up at their 
minimum. The different subplots show different contri- 
butions of the background (once again given by radia- 
tion). 
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FIG. 4: Regions where initial values for (0,o"i) will lead to a 
stabilisation of the field in the minimum of the potential, with 
a radiation background. The corresponding initial value for 
Sl r is indicated on each subplot. The plots are symmetric in 
o i so we only plot the positive values. The parameter values 
used in the potential are detailed in the text. 



FIG. 5: Stabilization regions for the initial conditions (</>, cr;) 
with a matter background. Initial value of Q m is shown for 
each subplot. The potential is the one used in figure 0] 



A few comments are in order. As can be seen from 
figures E| Hand El there is both a minimum and a max- 
imum value of <f> for which stabilisation at the minimum 
of the potential is successful. 

Let us begin with the left bound on the value of <f>. This 
comes from the type of potentials we are considering not 
being exactly an exponential potential. More precisely, 
they deviate most from an exponential shape for smaller 
values of oy. In these regions a scaling solution does not 
exist, and the evolution attractor is a scalar field domi- 
nated solution (the reasons for this are explained in more 
detail in Appendix 0). Consider a typical evolution like 
the one in figurcEl Initially the field will have values that 
correspond to a scalar field dominated (region 1) evolu- 
tion. Even if the background dominates initially (as in 
figure EJ, the field will eventually take over. In other 
words, there is a saturation point for the initial value of 
Pb, above which all evolutions will be similar (they only 
differ in a time shift, corresponding to the time during 
which the field remains frozen initially) . After some time 
the scalar field dominated solution becomes unstable and 
the field's energy density becomes kinetically dominated 
(region 2 in figure El- The background energy density 
will then dominate and the field will freeze at some con- 
stant value (region 3 in figure El- Note that the larger 
the initial energy density for the field (that is, the smaller 
the initial value of </>) , the longer it will take for the back- 
ground to freeze the field. For very small initial values, 
the field will run past the minimum before being frozen 
by the background. The limiting case corresponds to the 
field freezing to a value nearly coincident with the po- 
sition of the minimum (the case discussed in ref. [j]). 
A matter (7=1) type of background, having a slower 
redshift evolution, is much more effective at dominating 
the total energy density and freezing the field's evolution. 
This can clearly be seen in a much smaller left bound on 
the initial values of </> (see figure El- 




FIG. 6: Same results as in figure|lj but now fixing the initial 
value of p r instead of fixing Q r . The potential is the same. 
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The bound arising on the maximum value for (f> is also 
related to the scaling solution. Naively one would think 
that, as the initial value approaches the minimum, the 
field would acquire less kinetic energy and it would be 
easier to stabilize its evolution in the minimum. The 
reason we still get a bound on the maximum allowed 
initial value of (j>, is really due to the choice we made for 
the initial conditions. 

If the initial energy density for the field is higher than 
the scaling one, the field will have to reach the scaling so- 
lution in a kinetic dominated regime, and it will be easy 
for it to overshoot the minimum if it starts too close to it. 
As can be seen from figure for this type of potentials 
the scaling solution is essentially background dominated, 
so S7& is close to unity. By choosingto fix the value of 
f2b < 1 (which we did following ref. we immediately 
force the initial conditions for the scalar field to be away 
from the scaling solution. Furthermore, since the poten- 
tials are not exact exponentials, the scaling value actually 
increases and goes assymptotically to 1 as we approach 
the minimum. Therefore, for most initial values of fit,, 
the field will start above the scaling solution, leading to 
a right bound on the initial values of cf>. This bound will 
shift towards the minimum as we increase the initial fit,. 

Of course, eventually, as Qb approaches 1 we can get 
the right bound to coincide with the position of the min- 
imum. This is much easier to see if we choose to fix the 
initial value of pb, instead of f2&. In this case, the right 
bound always coincides with the position of the maxi- 
mum in the potential (just before the "roll over" point). 
The results are shown in figure El where the right bound 
is always the same for changing values of the initial pb- 

Given the dependence of the scalar potential on 
cos(acri) these plots are symmetric under the inversion 
Cj — > — &i (we only plot the positive Oi region). We can 
see that up to 85% of the parameter space in the imag- 
inary direction will lead to a stabilization of the field at 
the minimum. 



IV. KALLOSH-LINDE MODEL 

This model generalizes the original version of the 
KKLT model by admiting one additional componet in 
the superpotential. More explicitly, we write 

W = Wo + Aer aa + Be'' 3 ' 7 , (17) 

The particular example considered in ref. p"?| , sets pa- 
rameters, A = 1, B = -1.03, C = 0, a = 2tt/100, 
(3 = 27r/99 and Wq such that there is a supersymmet- 
ric minimum with zero cosmological constant, i.e. Wo 
is such that both W and F„ = K a W + W a vanish at 
some ay = cr C rit, &i = 0. Furthermore, there are a series 
of supersymmetric, AdS minima. All those are shown 
in figure [7| where we plot the scalar potential for this 
model, in the (ay, Cj) plane. The Minkowski solution is 
one of the two at = 0, and corresponds to the smallest 
value of oy. 




FIG. 7: Contour plot of the scalar potential, denned by 
eq. 117H . in the (oy, ay) plane. The values of the parameters 
are taken from ref. (17l | . see text. All minima are supersym- 
metric and AdS, apart from the one at <Ji — and smallest 
value of ay, which is Minkowski. 

Unfortunately, this choice of parameters makes it al- 
most impossible for the field to be stabilized if the initial 
value of the field is far from the minimum. In figure [S] we 
show the region of stabilization for a matter background, 
using a fixed p m for initial conditions. The reason for 
this is, again, related to the fact that the potential is 
very shallow for (f> < t/> m i n ~ 5.08, hence the scaling solu- 
tion is never the dominant attractor. 




FIG. 8: Stabilisation region for the initial conditions (4>,ai), 
with a matter background. The corresponding initial value of 
p m is shown on each subplot. We used the model of ref. |l7l[ . 
described in the text after eq. 1171 . 

Nonetheless, the mechanism of stabilising the field with 
a background fluid becomes impressive for models where 
the minumum lies within the range, in <j), of the scaling 
solution. One can see this by moving the minimum to 
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larger values of 0, for example, by decreasing B to B = 
-1.5. 

Following the approach of the previous section, in fig- 
ures [5] and ^| we show the bounds on the initial position 
of the fields that lead to successful stabilisation. As in 
the examples of sections IIII Al and IIII Bl the region of 
allowed initial conditions increases with the amount of 
initial energy density in the background. Also, as before, 
a smaller background equation of state will lead to better 
stabilization. 

In the imaginary direction, this type of potential is 
much more difficult to stabilize than the ones used in the 
previous section. For these examples, we can only cover 
around 0.6% of the available range of initial conditions. 
The range of values for at are similar in both cases, so this 
is only due to the fact that the examples in this section 
have a much larger period in the imaginary direction. 
This can be easily checked by comparing figures \7\ and ^ 
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FIG. 9: Region of the initial values in ((f>, Oi) for which stabil- 
isation at the minimum occurs in the model of ref. |l7|. with 
parameter B set to —1.5. We used a radiation background, 
fixing the initial value of p r . 



a = 4.0 xlO" 10 k; 4 



PmO = 1.0 xlO^ Kp 
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FIG. 10: Region of the initial values in (<f),<Ji) for which 
stabilisation at the minimum occurs in the model of ref. [TUl . 
with parameter B set to —1.5. We used a matter background, 
fixing the initial value of p m . 



the moduli fields. We have evaluated the region of al- 
lowed initial conditions for the fields and confirmed that 
it increases for larger values of the initial background en- 
ergy density po and lower values of its equation of state. 
The effectiveness of this mechanism is also dependent 
on the existence of a scaling solution as the attractor. 
Regions in which the potential is considerably shallower 
than an exponential lead to a scalar field dominated at- 
tractor, working against successful stabilisation in the 
minimum. 

In the second part of this paper we have turned our 
attention to analyse the Kallosh-Linde model of ref. ^3] • 
The conclusions are qualitatively equivalent but we note 
that only a small interval of initial conditions are allowed 
in the imaginary direction, making it more difficult to 
obtain stable moduli solutions in this case. 



APPENDIX A 



V. CONCLUSIONS 

In this article we have studied moduli evolution in the 
context of the so-called KKLT-like models, that arise out 
of Type IIB string theory, where fluxes and the effect of 
anti-D3 or D7 branes have been taken into account. We 
have mainly looked at the allowed region of initial con- 
ditions of the moduli fields which leads to their eventual 
stabilisation at the supersymmetric minima of the scalar 
potential. Employing an approach first developed in Q, 
we found that a background perfect fluid has the ability 
to slow down the fields preventing them from running 
past the minimum, in agreement with Q. 

Moreover, we have extended the work of ref. £| to in- 
clude the dynamical evolution of the imaginary part of 



In this appendix we summarise the key results pre- 
sented in and [2l|, and describe how they can be ap- 
plied to the discussion in the main text. 

In the first of these papers , we considered the issue 
of stabilising the dilaton in the context of supcrstring 
cosmology without the explicit use of non-trivial three 
form fluxes, but including gaugino condensates and ex- 
ternal sources such as radiation. Making use of the fact 
that in these models the scalar potentials arising out of 
gaugino condensates are generally exponential in nature, 
we demonstrated how sca ling solutions associated with 
exponential potentials [Islll9l,l20j could be slightly mod- 
ified for these models, but basically had similar outcomes. 
In particular we showed how, due to the friction of the 
background expansion, the energy density in the scalar 
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field could become, for a period, a fixed fraction of the 
background density , allowing the field to be trapped 
in the minimum of its potential, generally as it left its 
scaling regime. The key point was that previous scaling 
solutions available for pure exponential potentials were 
also approximately valid for quasi-exponential potentials. 
This was further reinforced in |2l| . where a phase space 
analysis of generic quasi-exponential potentials was car- 
ried out, looking for scaling solutions. 

Consider a scalar field <fi, with canonical kinetic terms, 
evolving in a FRW Universe containing a fluid with 
barotropic equation of state p 7 = (7 — It is use- 

ful to analyse the system in terms of the new variables 
x = J^ H and y = ^7. , as they allow for the determi- 
nation of the scaling solutions. The evolution equations 
can then be written as 12(1 (23 



(Al) 



H' = -°-H[2x 2 +1 {l-x 2 ~y 2 )] 
x' = -3x + \^y 2 + ^x[2x 2 + 1 (l~x 2 -y 2 J\A2) 

V = -\\J\xy+ i -y[2x 2 + 1 {l-x 2 ~y 2 )] (A3) 

(A4) 



A' = -V6A 2 (r-l); 



where a prime denotes a derivative with respect to N 



In — - — , and we defined two new variables 

^initial ' 



pect to 

MM, 



1 dV 

KpV dtp 



r - 1 



KpX 



(A5) 



For the pure exponential potential, A is a constant and 
r = 1, so in a sense T is measuring how much our poten- 
tial deviates from an exponential one. These equations 
can be solved easily in a number of regimes 0, |21| . The 
types of regimes available were already described in sec- 
tion and are shown in figure [21 Here we will present 
a more analytic description. 

Let us start with the kinetic dominated regime, that is 
region 2 in figure In this region, the evolution is dom- 
inated by the fields kinetic energy and so the potential 
energy, y, can be neglected in the equations of motion 
eqs. l|A2HA3j) . This leads to the solution 



x = 1 



1 



T 2 

x 3(2- 



7 )JV 



-1/2 



(A6) 



where xq is the initial condition for x (taken at N = 
0, for simplicity). Notice that, since we are neglecting 
the potential terms (y), this solution is really potential 
independent and is valid as long as we are in a kinetic 
dominated regime. Clearly the kinetic energy of the field 
which is proportional to x 2 is decreasing rapidly with 
time. The solution for <f> follows and is given by 



sinh 



•''0 



-3(2- 7 )AT/2 



, (A7) 



where 4>o is the initial value of the field. At early times 
the solution corresponds precisely to that of eq. (8) of Q, 
where the energy density is dominated by the kinetic en- 
ergy of the scalar field. As time carries on though, and 
the field slows down, it eventually tends to a constant 
solution as can be seen in eq. l|A7p . This constant value 
corresponds to region 3 of figure [3 

Let us now turn our attention to the scaling evolution, 
that is region 4 in figure |3 In [2lJ an approximate equa- 
tion of state for the field evolution in the scaling regime 
was derived. This is given by 



70 



7 +(2r-l)A 2 /3] 
y^+( 2r - 1 ) A2 /3] 2 + 87(r-l)AV3. 



(A8) 



When r — > 1 (the exponential case), 7^ tends to 7, that 
is one recovers the usual result that the field evolution 
mimics the background. Moreover, this scaling solution 
is only stable if it satisfies the following conditions, 



7* < 
14> < 



7A 2 

~T' 

3 7 - 2A 2 
6 



and 



12A 2 (2 + 7 ) 
(3 7 -2A 2 ) 2 



(A9) 



In the particular case of double exponential potentials, 
V(4>) = Vo exp(ae & ^), it is possible to obtain an approx- 
imate solution for the field evolution. The most useful 
way to present it is in a recursive form 0, l2lj 



S 
Si 
S k 



3 Po 7(2 - 7) 
2 V a 2 b 2 



lln 

a 



a 

S --]n(Sk-i)-—N, 



(A10) 



where S = exp(60). It is easy to check that in the regions 
where the potentials discussed in sections IlIII and llVl are 
close to an exponential in ay, these solutions are a very 
good fit to the scaling evolution with S ~ oy . 

Finally, there is a third possible type of evolution, 
shown in region 1 of figure |2 the scalar field dominated 
solution. Its effective equation of state was shown in [2l| 
to be 



14> = ^*/3 : 



(All) 



where we have defined 



r(N) =0o + 



2^6 

3(2-7) 



sinh 1 



•' : o 



1- y/1 - 4(r - l)A 2 /3 
0^1)A 



(A12) 
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for r^l, and A^ = A otherwise. Again, a condition for ACKNOWLEDGMENTS 
the stability of the solution was derived, and is given by 

6 (A -A)-^ + i(A^ + AA -6 7 ) < 0, (A13) 

3 - 6 A + ^A43A - 2A) < 0. (A14) 

<\tj> z 

We can now apply this knowledge to the potentials we 
are dealing with in this paper. In figure ITU we compare 
the evolution of the equation of state resulting from a nu- 
merical integration with the scaling and scalar field dom- 
inated solutions given above, in their region of stability. 
The model used is the KKLT example presented in sec- 
tion [ffij also used to produce figure |2 We can clearly see 
from figure ITT1 that in this specific example the attractor 
solution only makes sense for </> > 3.7. Therefore, if the 
field starts its evolution before this value, it is either in 
the scalar field dominated solution or in a kinetic energy 
dominated solution. This makes the stabilisation of the 
field more difficult to achieve for smaller initial values of 
the field. 

2, 1 , 1 , r-^ r-, 1 



1.5- 




0.5 



° -4 -2 2 4 6 

FIG. 11: Dependence of the scalar field equation of state 
with the value of the field. Solid line: the numerical integra- 
tion; dash-dot line: scalar field dominated solution; dash line: 
scaling solution. The model is the same we used in section 
IIHI with a background of dust (7 = 1), i.e. the same as in 
figure 2. 
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